Unique Fock quantization of scalar cosmological perturbations 
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We investigate the ambiguities in the Fock quantization of the scalar perturbations of a Friedmann- 
Lemaitre- Robertson- Walker model with a massive scalar field as matter content. We consider the 
case of compact spatial sections (thus avoiding infrared divergences), with the topology of a three- 
sphere. After expanding the perturbations in series of eigenfunctions of the Laplace-Beltrami op- 

Jvj ' erator, the Hamiltonian of the system is written up to quadratic order in them. We fix the gauge 

of the local degrees of freedom in two different ways, reaching in both cases the same qualitative 
results. A canonical transformation, which includes the scaling of the matter field perturbations 

p^i ' by the scale factor of the geometry, is performed in order to arrive at a convenient formulation 

of the system. We then study the quantization of these perturbations in the classical background 
determined by the homogeneous variables. Based on previous work, we introduce a Fock represen- 
tation for the perturbations in which: (a) the complex structure is invariant under the isometries 
of the spatial sections and (b) the field dynamics is implemented as a unitary operator. These two 
properties select not only a unique unitary equivalence class of representations, but also a preferred 
field description, picking up a canonical pair of field variables among all those that can be obtained 
by means of a time-dependent scaling of the matter field (completed into a linear canonical trans- 

' r\ ' ■ formation) . Finally, we present an equivalent quantization constructed in terms of gauge-invariant 

pjH 1 quantities. We prove that this quantization can be attained by a mode-by-mode time-dependent 

linear canonical transformation which admits a unitary implementation, so that it is also uniquely 
determined. 
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I. INTRODUCTION 

The theory of cosmological perturbations plays a major role in the explanation of structure formation in the uni- 
verse, describing the growth of originally small perturbations in an otherwise homogeneous and isotropic background 
up to the nonlinear regime. Lifschitz [lj was the first to study the perturbations of the metric tensor in a Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model. However, it was soon realized that the interpretation of this kind of 
perturbations was obscured by the freedom to perform gauge transformations, i.e., infinitesimal coordinate trans- 
formations that keep fixed the background spacetime. Two different ways of avoiding this problem were proposed. 
On the one hand, Hawking Q developed the first attempts to construct a fully covariant, gauge-independent theory, 
which was completed by Olson for the case of an isentropic perfect fluid with vanishing spatial curvature [3[ . These 
works served as starting point for the covariant formulation of Ellis and Bruni [4,] . On the other hand, Bardeen [5( 
showed how to combine the metric and matter perturbations so as to obtain gauge-invariant quantities. This latter 
approach has been extensively studied [6( . 

Despite the fact that perturbation theory in cosmology accounts for the growth of the inhomogeneities that should 
eventually produce the large-scale structure of the universe, by itself it does not clarify the origin of those perturba- 
tions. At this point, additional hypotheses must be accepted. The paradigm of cosmological inflation [7j provides a 
mechanism for generating inhomogeneities and, moreover, solves certain problems of the standard model of hot big- 
bang cosmology, such as the horizon and flatness problems [8[ . According to inflation theory, the universe underwent 
a period of very rapid expansion in its early history, which is generally attributed to the presence of a hypothetical 
scalar field: the inflaton. The quantum fluctuations of this field give rise to the primordial perturbations. Among 
the many models of inflation that have been proposed, one of the most popular is chaotic inflation, put forward by 
Lindc 9]. In this model, the inflaton, subject to a potential — which can be as simple as a quadratic one, owing to 
its mass — , is chaotically distributed before inflation. The patches of the universe where the field has large values 
expand exponentially, becoming nearly homogeneous in the process. One of these regions would encompass the entire 
observable universe. 

In order to satisfactorily understand structure formation in the universe, therefore, a quantum or at least semiclas- 
sical description of inflation is needed. In the absence of a complete quantum theory of gravity, one is forced to adopt 
the latter approach and try to quantize the inhomogeneities in a classical, homogeneous and isotropic background. 
However, this approximation suffers from a series of problems, like for instance the generic ambiguities that plague 
Quantum Field Theory in Curved Spacetime (QFTCS). These ambiguities appear even in the case of linear fields, 
for which the standard procedure of Fock quantization is available. To begin with, it is far from clear which are the 
variables that one should regard as basic ones. Indeed, inasmuch as the dynamics of the background is treated in 
a totally different way compared to those of the field which is being quantized, a splitting between them must be 
defined. For instance, a scaling of the variables by a time-dependent background function (such as the scale factor 
in an FLRW universe) certainly modifies their dynamics, while preserving the linearity of the system. When there 
is a finite number of degrees of freedom, the classical transformations of this type (linear but time dependent) are 
generally implemented as unitary transformations in the quantum theory. Hence, the quantization of the original and 
of the scaled variables are unitarily equivalent, and thus they lead to the same physical predictions. By contrast, 
that is not necessarily the case in a field theory. Moreover, there is another, more standard source of ambiguity 
in the quantization of linear fields, inherent to QFTCS. It is the election of a Fock representation of the Canonical 
Commutation Relations (CCR's). Once again, this problem does not appear in systems with finite degrees of freedom, 
where the Stone-von Neumann theorem [lOj guarantees the uniqueness (up to unitary equivalence) of the (strongly 
continuous, irreducible, and unitary) representation of the CCR's in their Weyl form. But, when dealing with an 
infinite number of degrees of freedom, infinitely many inequivalent representations are possible. This ambiguity is en- 
coded in the choice of a complex structure on the phase space T [ll|, [l2| • A complex structure J is a linear map whose 
square equals minus the identity (i.e., J 2 = —I). We also require its compatibility with the symplectic structure fi, in 
the sense that it must be a symplectomorphism and the bilinear map f2( J-, •) must be positive-definite. In this way, 
[H(J-, •) — iH(-, -)]/2 provides an inner product, which endows T with the structure of a complex Hilbert space — the 
one-particle Hilbert space. Equivalently, one can see that a complex structure J decomposes the complexification of F 
into two eigenspaces of eigenvalues ±i, thus generalizing the concepts of positive and negative frequency available in 
stationary scenarios. With all these ingredients, it is possible to define the full Hilbert space of the theory, a vacuum 
state in it, and finally, a representation of the CCR's. 

In spite of these ambiguities, there exist criteria to select a preferred quantization in some specific situations. A 
clear example is Minkowski spacetime, where Poincare symmetry permits one to select a complex structure for a given 
Klein-Gordon field. In systems with a lower degree of symmetry, other kind of conditions are required, such as energy 
considerations in the case of stationary spacetimes, where a time- like Killing vector is available [13}, or the existence of 
an invariant Gaussian solution to a regulated Schrddinger equation in the specific case of a de Sitter background [14| . 
Furthermore, in the absence of stationarity, it has been proven recently that the unitary implementation of the field 



dynamics is a useful criterion to select a unique Fock representation of the CCR's (combined with the invariance 
of the vacuum under any spatial symmetry remaining in the field equations). For instance, in the Gowdy models 
(inhomogeneous cosmological models with compact spatial sections and two space- like Killing vectors [HI]), after a 
symmetry reduction of the system, the inhomogeneities can be treated in a formally identical manner as a scalar 
field with time-dependent mass in an auxiliary spacetime of lower dimension. The spatial sections of the auxiliary 
spacetime can be circles or two-spheres, depending on the topology of the spatial sections of the Gowdy model into 
consideration. In both cases, the complex structure that would be natural if the field were massless allows for a 
unitary quantum implementation of the field dynamics and, besides, shares the symmetry of the spatial sections 
of the auxiliary spacetime (which the vacuum state inherits). In fact, these two properties select a unique uni tary 
equivalence class of Fock representations: all the representations with these attributes are unitarily equivalent jlq - fl8J . 
Furthermore, if we change the field description by performing a linear canonical transformation that scales the field by 
a time-dependent factor, unitarity is lost: there is no Fock representation for the new canonical pair of field variables 
in which dynamics is represented in terms of a unitary operator [19| . These results also apply in the case of a scalar 
field with time-dependent mass on the three-sphere [20, [21| , that is to say, in a spacetime whose spatial sections are 
three-dimensional spheres, a situation which is especially relevant for our discussion in this article. 

The aim of the present work is to discuss in detail the application of the requirements of unitary evolution and 
symmetric vacuum state in order to remove the ambiguities in the Fock quantization of scalar fields in a context of 
particular physical interest: the scalar perturbations in an FLRW model with a massive scalar as matter content and 
spatial sections with the topology of a three-sphere. In Sec. UH we analyze the system classically. Its Hamiltonian 
is derived in Sec. Ill Al briefly reviewing the work of Halliwell and Hawking [22|. In doing this, we expand the 
perturbations in the modes of the Laplace-Beltrami operator and conduct a perturbative analysis on them. Then, in 
Sec. Ill Bt we fix the gauge of the local degrees of freedom and perform a canonical transformation of the remaining 
variables so as to reach an adequate formulation of the system, especially convenient for the rest of the analysis. 
In particular, we scale the perturbation of the matter field by the FLRW scale factor. We conclude the classical 
treatment writing the dynamical equations of the new canonical pair of field variables. 

We study the Fock quantization of the scalar field perturbation in the classical background determined by the 
homogeneous variables in Sec. Mil In Sec. 1111 Al we introduce a Fock representation that fulfills the requirements 
of symmetry invariance and unitary dynamics. Any other representation with the same properties is shown to 
be equivalent to it in Sec. IIIIB1 The rest of Sec. Ill is devoted to the analysis of different descriptions of the 
system, with other elections of canonical variables. First, in Sec. MI CI we consider a time-dependent scaling of the 
perturbation as a whole, completed into a linear canonical transformation, and prove that the new dynamics cannot be 
implemented unitarily. More general mode-by-mode transformations are then taken into account in Sec. IIIID1 where 
the system is reformulated so that a natural quantization with unitary evolution is possible. This new quantization 
is unitarily equivalent to the one previously introduced, and the Hamiltonian is such that no cross term between the 
field configuration and the field momentum appears. 

In Sec. IIV1 we define a set of gauge-invariant quantities and write their expressions in the chosen gauge. The 
dynamics of these gauge invariants is seen to be equivalent to that of a Klein-Gordon field with time-dependent mass. 
We show that any Fock representation which admits a unitary implementation of this dynamics is unitarily equivalent 
to the one we already had. Finally, a different gauge choice is explored in the Appendix. 



II. THE CLASSICAL MODEL 
A. Hamiltonian of the system 

We assume a globally hyperbolic spacetime and start with the ADM form of the metric [23[ . We use the standard 
notation N for the lapse function, N % for the shift vector, and hij for the three-metric on the spatial sections of 
constant time t. Spatial indices are denoted with lowercase Latin letters from the middle of the alphabet, and are 
lowered and raised with the metric hij and its inverse, respectively. This decomposition allows one to express the 
Einstcin-Hilbert action for General Relativity as the time integral of the Lagrangian [2J] 

L g = ^ Jd 3 x NVh [KijK* - (K\) 2 + 3 R] , (1) 

where h is the determinant of the three-metric, K%j is the extrinsic curvature of the spatial sections, 3 R is their scalar 
curvature, and G is the Newton constant. 

The total Lagrangian of the system is the sum of L g and the contribution of the matter content: a scalar field $ of 



mass to. The latter can be expressed in terms of the ADM variables as [22J 

^) -2N i d i ^^--(N 2 h ij -^N^d^dj^-N 2 ^ 2 ^ 2 . (2) 



i r~ Vh 



L m -2 dx w 



where di is the flat derivative with respect to the spatial coordinate x % . 

In the FLRW model, N depends only on time and N l vanishes, whereas the spatial sections have constant curvature. 
This is no longer the case if inhomogeneities are present. The description of such inhomogeneities in an FLRW 
background is greatly simplified by the introduction of a convenient basis of modes on the spatial sections. For the 
three-sphere S 3 , a natural basis of scalar functions (complete in the space of square-integrable functions with respect 
to the volume element provided by the three-metric) is given by the eigenfunctions of the Laplace-Beltrami operator, 
namely, the hyperspherical harmonics [22, 25, 26]. The corresponding eigenvalue equation can be written in the form 

fi^ViVjQPm = -<4QL, (3) 

where the symbol V indicates covariant derivatives with respect to the round metric on S 3 , which we denote by Clij. 
Note that its inverse is represented here by Q y , and that in FLRW we have hij — a 2 (t)ilij, with a(t) being the scale 
factor of the spatial geometry. 

The eigenvalues satisfy the relation w„ = n — 1, where n £ N + can be any positive integer. The corresponding 
eigenspaces have a degeneracy equal to n 2 , since the labels I and m can take all integer values in the ranges < I < n—1 
and —l<m<l. For convenience, we adopt a basis of real hyperspherical harmonics, because we are dealing with 
a real scalar field (if one started with complex harmonics, real ones can be straightforwardly obtained by taking the 
real and imaginary parts of the basis functions). Moreover, the hyperspherical harmonics are normalized so that they 
are orthonormal, in the sense that 

fx 



ttQ? m Q?> m > = <5™" S u >S mm >. (4) 

In the following, the labels n, I, and m will be collectively denoted by n, except if it is more convenient to show them 
explicitly. 

Vector and tensor harmonics can be obtained from the scalar ones by means of covariant derivatives: 



3" 



(Pi) n = — ViQ n , (P t3 ) n =V J (P i r + -n iJ Q n , ng{2,3,...}. (5a) 



These harmonics are not complete in the respective spaces of vector and tensor fields on S 3 , as there exist other vector 
and tensor harmonics which cannot be obtained in this way. Nevertheless, we will consider only inhomogeneities 
constructed from scalar perturbations (and their derivatives) . A complete analysis shows that the dynamics of such 
perturbations is completely decoupled from the genuinely vectorial and tensorial ones [22j . 

The variables of the model can be expanded in series of hyperspherical harmonics, the n = 1 term corresponding 
to the homogeneous variables of the FLRW model. We can make this explicit by writing 



(G) 



hij =cr 2 e 2Q ^- f l + 2V2ireJ2 a ™Q n ) + cr 2 e 2a 6V2iTe ^ b n (Pij) n , 

\ n / n 

N = aN U + V2TreJ29 n Q n ), N t = a 2 e«v^6^ *„(£)", $ = \ ( "4-P + e E /"0" J ' 

\ n / n \ * n / 

All these sums start in n — 2, except the metric contribution with coefficients b n , which begins in n — 3 because 
the corresponding tensor harmonic with n = 2 turns out to be proportional to flij. The factor a = y2G/37T has 
been introduced for convenience. The "background" functions Nq, a, and ip parameterize the homogeneous solutions, 
whereas the coefficients a n , b n , g n , k n , and f n describe the inhomogeneities. They are all functions of time only, since 
the hyperspherical harmonics absorb all the spatial dependence (apart from that of the components of the metric f2jj ) . 
We assume that the inhomogeneities are sufficiently small so as to justify a perturbative analysis up to second order in 
the parameter e (namely, up to second order in the inhomogeneities). Note that if all the perturbations are neglected, 
the FLRW model is recovered, with crJVo being the homogeneous lapse function, ae a the scale factor, and {-\/2'Ka)^ 1 ip 
being a homogeneous scalar field of mass m = raj a. 



Introducing expressions (0 in Eqs. (QJ and ([2]), one can obtain the Lagrangian of the system within the perturbative 
framework [22f . A Legendre transform with respect to the time derivatives of the variables a, ip, a n , b n , and f n leads 
to a Hamiltonian of the form 



H = N 



tflo + ^El^+^l" 



e 2 yk n H^ + 0(e 3 ), 






(7) 



which is a linear combination of constraints, with Nq, Nog n , and k n acting as Lagrange multipliers. The subscripts 
of the quantities appearing in the last expression indicate their association with the lapse function (|) or the shift 
vector (_), as well as their perturbative order in terms of the inhomogeneous modes (namely, their degree as polynomials 
of the coefficients of the inhomogeneities) . 

Calling generically ir q the momentum conjugate to any variable q [27| . the zeroth-order Hamiltonian is written as 



#|o = \e- 3a (-irl + tt 2 + e 6a m 2 ip 2 - e Aa ) , 



while the higher-order terms have the explicit form 



(8) 
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n l - 1 
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„4a J 1 /„2 5\ „2 , 1 L2 



™ 2 -4, 2| 2 /2 



(» 2 -l)«n+i(" 2 - 7 )^n + ^(» 2 -4)^n-(n 2 -l)^ 



< '"'nr | ^a 2 , - 6^— ^ 2 £ + 6y>a„/ B + /* 



11 2 



i?" = -e- a 
- 1 3 



»t£ - 3a n irl - 2TT a Tr an + 27t v 7t /ti - -e ia [(n 2 + ±) a„ + (n 2 - 4) 6„] 



e ba m^ (3y>a„ + 2/„) 



n 2 - 4 
n — 1 



ft ™ + 4 l2 i"" n ) ^ + SfnKtp - 7T Qti + TT K 



(9) 



(10) 
(11) 



B. Gauge fixing and reformulation of the system 



The presence of constraints indicates that the system includes unphysical degrees of freedom. Here, we partially 
fix this gauge freedom by adopting the longitudinal gauge. A different gauge choice is addressed in the Appendix, to 
show that the results are not sensitive to the particular gauge adopted. 

More specifically, we want to remove the freedom associated with the Lagrange multipliers g n and k n (corresponding 
to the constraints H," = and H™ = 0, respectively), setting Ni = and h%j ex Oy. To this end, we impose 



0, C n 



^rh^OL dJnTTu 



0. 



(12) 



These are acceptable gauge- fixing conditions as far as n > 2, since in this case their Poisson brackets with H n x and 
if,™ satisfy 



-i,tY{J-Sl {^'SlW(n 2 -4)>0. 



(13) 



The symbol ~ indicates that the identity holds when one makes use of the constraints and the gauge-fixing conditions. 



Eqs. (|12l) . together with the constraints, allow one to write a n and nb n in terms of /„, Wf n , and the homogeneous 
variables. In particular, 

whereas 7T(, n = 0. Besides, consistency of the gauge-fixing conditions with the dynamical evolution requires the 
imposition of the following additional restrictions: 

= {b n ,H}n~ e - a k n + O(e), = {C n , H} w i(n 2 - 4)N e a (a n + g n ) + 0(e). (15) 

Therefore, we get k n — and g n = —a n up to higher-order perturbative corrections. 

In the case n = 2, conditions (fl~2]) are no longer applicable since, by construction, the coefficients b n exist only for 
n > 2. Actually, for the second mode, one can always fix the gauge so that a-iim — film — 0. This leaves no remnant 
physical degree of freedom for this mode. 

After this gauge fixation, the spacetime metric reads 



ds 2 = -a 2 N 2 ( 1 - 2V27reJ2 a nQ n ) dt 2 + o 2 e 2& Sl l3 ( 1 + 2\/2 7re^ a n Q n ) dx i dx j + 0(e 2 ). 



(16) 



In the reduction of the system, it is convenient to introduce new coordinates on phase space so as to reach a 
canonical set with respect to the new, reduced symplectic structure. Namely, we keep 7r Q , n v , and /„, while their 
conjugate variables [27| become 

a = a + —^a 2 l , (p = ip + 3e 2 ^ a nfn, tt/„ = tt/„ - 3a„7r v . (17) 

n n 

The reduced Hamiltonian can be written as 

H = N (H l0 + t 2 Y.^)+O(e 3 ). (18) 

Since the new homogeneous variables differ from the old ones only by quadratic terms in the perturbations, the 
expression of the zeroth-order Hamiltonian is unaffected by the transformation: it is still given by the right-hand side 
of Eq. ([5]). provided that the old variables are replaced with the new ones. As for the second-order Hamiltonian, it is 
of the form 

#|2 = \^ (KA + 2E]J n n fn + E] f f 2 ) , (19) 

where we have introduced the notation 

E^ = e- 2& (l-^-nl), (20a) 



n 2 -4 
4~ 



% = ^rr7 e_6 "^ ( 37r «^ - e 65 m 2 ^) , (20b) 

E] f = e 2& (n 2 - 1) + e 4 "m 2 - 9e- 25 7r 2 - -^— e - && (3^ Q 7r v - e 6 "m 2 ^) 2 . (20c) 

Before concluding the classical treatment of the system, we perform a canonical transformation that reformulates 
the system in a especially suitable form. First of all, we scale the perturbation of the matter field by the FLRW scale 
factor: /„ h-> e a f n . This transformation (that we will call the Mukhanov scaling; see, e.g., Ref. [28]) is frequently 
adopted when working in an FLRW background, as it simplifies the equations of motion. We can complete this change 
into a canonical transformation by introducing the inverse scaling in the conjugate field momentum. Besides, we allow 
also for a contribution to this new momentum that is linear in the field configuration, with the aim at removing cross 
terms in the Hamiltonian which would couple the configuration and momentum variables of the perturbations — at 
least up to subdominant terms in the large-n limit. Apart from further simplifying the system, in this way we adapt its 



description to the most convenient form for the asymptotic analysis of the forthcoming section. These new variables, 
which have canonical Poisson brackets at the considered order in the perturbations, are 

5 = 5+ ?E^' ^a = 7T Q + e 2 ^ (-/„tt/„ + tt q /^) , f n = e a f n , Kf n = e~ a (TTf n - TT a f n ) , (21) 

n n 

while ip and ir^ are left unchanged. 

In terms of these new variables, the Hamiltonian has the same structure as in Eqs. (|18|) and (|19j) . with the functional 
form of the zeroth-order Hamiltonian still being given by Eq. ©. On the other hand, the coefficients of the second- 
order Hamiltonian are now 

K, = 1 - -^re- 4a 4, (22a) 



'<p> 
3 

n 2 - 4 ' 



E% = — je- 6a 7r v (27f a 7r v - e 6a m 2 ^) , (22b) 

% - n 2 - ± + ±e aa m a (2 - 3^ 2 ) - ^e~ 4a (tt 2 + 15tt 2 ) - — — e~ 8a (2^ - e 6a m 2 ^) z . (22c) 



Thus, as we had anticipated, we have implemented the Mukhanov scaling while retaining the behavior of the cross 
terms in the Hamiltonian, ij? = 0(n~ 2 ), again like E 1 }^ in Eq. (|20bl) . 

On the other hand, the expression of a n as a function of these variables is 

a n = -^— [e~ 5a (e 6s m 2 p - 2n a n v ) f n + e- 3 %,7f /n ] . (23) 



Taking this into account, the spacetimc metric can be easily obtained via Eq. (JT€ 

To end this section, we will write the dynamical equations of the system. In doing this, we will use the conformal 
time 77, defined by N dt — e a dr\. Derivatives with respect to this time will be denoted by a dot. Hamilton's equations 
give 



e 



-2o 



a = -e~ 2a Ti a + 0{e 2 ), tt & = — ( - 3^ 2 + 3tt£ - 3e ba mV + e 4a ) + 0(e 2 ), (24a) 

$ = e- 2a TT v + 0{e 2 ), n v = -e 4a m 2 p + 0(e 2 ), (24b) 

/„ = {K^ fn + %/„) + 0(e) , (24c) 

# /n = - (%7f /n + %/„) + 0(6). (24d) 

Using Eq. (|24cl) . we can isolate 7fj n and write it as 

7T/„ = (1 +Pn)fn + Qnfn, (25) 

where we have neglected higher-order contributions and introduced the time-dependent coefficients 

v - 3 ^ „ - -3e- 2& n 2 ^ y - e 6a m 2 ^ 

P " ~ (n 2 - 4)e 4a - 3^ 2 ' 9 " ~ v (n 2 - 4)e 4a - 3tt 2 ' l ' 

Combining the two last equations in the set (I24p , we obtain the equation of motion of f n : 

J^E}\ ~ E% ~ (E%) + E] f E^ f n - ^f n + f n = 0(e), (27) 

7T7T / 7T7T 

which can be written in the simpler form (neglecting again higher-order terms) 

/»+ r„/„ + (w* +«„)/„= 0, (28) 



where 



27r a 7r ¥ , — e 6a m 2 <f 
(n 2 - 4)e 4a - 3tt 2 



6e~ 2a 7r V - — y (29al 

1 ,* , 1 



•s„ = - e 2a rn 2 - ± e - 4a (tt 2 + 21tt 2 



- ^e _4S (*? + 21^ 2 + 3e 6a m^ 2 ) + O (n~ 2 ) . (29b) 



III. FOCK QUANTIZATION OF THE PERTURBATIONS 

Let us study now the Fock quantization of the scalar field perturbation in the background determined by the 
homogeneous variables, which we will describe classically (to all purposes, we will treat them as mere functions of time). 
This semiclassical treatment should provide a good approximation to the behavior of the system as long as quantum 
effects on the background do not become important and one can disregard non-conventional quantum geometry 
phenomena (e.g., those of polymeric type predicated by Loop Quantum Gravity |29l |30() in the inhomogeneities. 
Thus, we focus on the (standard) quantum fluctuations of only the local degrees of freedom, neglecting quantum 
corrections that could arise from the zero modes parametrized by the variables a, <p, and their momenta. 

The arguments presented in the following two subsections are similar to those detailed in Ref. [2(| for the Fock 
quantization of a Klein-Gordon field with a time-dependent mass on S 3 . Despite the differences between the two 
situations, the analysis can be adapted without major difficulties. In order to make this manuscript self contained, 
we now summarize the master lines of the proof discussed in Ref. [20( , and explain the modifications needed for our 
specific model. 

A. A quantization with unitary dynamics 

As we have commented in the Introduction, QFTCS is plagued with ambiguities. These include the choice of a 
field description. Among all those that can be reached by a time-dependent scaling of the field, we will initially select 
the canonical pair formed by the perturbative modes f n and their conjugate momenta %t , although we will discuss 
other possible options later on. In the previous section, we derived the equation of motion (J28I) for f n . In particular, 
note that all the modes with the same value of n evolve in the same way. Besides, in the asymptotic limit n — >• oo, 
one recovers the equation for the modes of a free massless scalar field on S 3 , at least apparently, because the function 
r n is negligible [see Eq. (|29a|) ]. and the function s n [in Eq. (|29b[) ] is much smaller than w^. On the other hand, the 
momentum Wt obeys Eq. (|2"5)l . and hence differs from the time derivative of /„ only by terms of the order of -nT 2 . 

Even after selecting a canonical pair, there is still an infinite ambiguity in the choice of a Fock representation for the 
corresponding CCR's. This freedom in the construction of a Fock representation amounts to the selection of a complex 
structure that is compatible with the symplectic structure. To fix this complex structure, let us first introduce the 
annihilation and creation variables that would be natural to adopt in the case of a free massless field: 



*/„ \ 1 (u n i \ ( f„ 



a/, I ^/27~ ^ 7 V -'„ - •' 



(30) 



with the symbol * standing for complex conjugation. Then, we choose the complex structure Jo that is diagonal in 
this basis of phase space variables: 

*C£)-(^)C£)- 

In the following, we will call jo the 2x2 matrix that appears in this last equation. 

The complex structure introduced above shares the symmetry of the spatial sections, i.e., it is invariant under SO (4) 
rotations. It is clear that such transformations do not mix variables with different label n, as they commute with 
the Laplace-Beltrami operator. Moreover, the configuration and momentum spaces associated with a given fixed n, 
namely the sets 

J"„ = {/„i m ; fixed n}, V n = {nj nlm ; fixed n}, (32) 

transform with the same irreducible representation of SO (4) [20(, so creation variables do not mix with annihilation 
variables under the kind of transformations that we are considering. 

Owing to the fact that the modes are dynamically decoupled, the symplectomorphism IA that gives their time 
evolution is block diagonal. Thus, fixing a(n arbitrary) reference time 770, we can write 

af>)j= w »(^°)yj> )j- (33) 

The matrix lA n depends on n, but not on the degeneracy labels I and m, because the evolution is insensitive to them. 
Besides, since U n represents a real transformation, it has the general form 

U n {r,,m)=( a R &^\ ^MV (34) 



Furthermore, as the dynamics is a symplectomorphism, the Bogoliubov coefficients a n and j3 n satisfy the relation 
|a n | 2 — \P n \ 2 — 1, so that Poisson brackets remain invariant under the transformation. 

In the quantum theory, annihilation and creation variables a n and a n are promoted to operators a n and a n , whose 
evolution is given by an analogue of Eq. (1331) : 

M-^^'W')- (35) 

The complex structure Jo has the particular property that the Hilbert-space operator that implements the above 
evolution exists and is unitary, as we will now prove. The necessary and sufficient condition for that to happen is that 

J -U(r),r}o)JoU- l (r),r)o) (36) 

must be a Hilbert-Schmidt operator (for all times r}) on the one-particle Hilbert space defined by Jo [3l|, 133 . This 
condition guarantees that the evolution of the vacuum state at 770 is a (normalizable) state at all times, and can 
be rephrased as the requirement that the beta-coefficients of the dynamical transformation are square summable (in 
complex norm), namely: 

^|^ n (r ? ,j 7o )| 2 <oo, (37) 

n 

which, recalling the n 2 -fold degeneracy of each eigenvalue n, is equivalent to demanding the square summability of the 
sequence {n/3 n } n . Consequently, the unitary implementation of the dynamics relies only on the asymptotic behavior 
of the Bogoliubov coefficients j3 n for large n, which can be deduced from the equation of motion (|28[) . In order to 
study the asymptotic limit of the solutions, we write them in the form 

f n (v) = A n e^ e ^) + A* n e^ e '^\ (38) 

where A n is a complex parameter that is determined by the initial conditions on f n (i.e., the initial values of the 
mode and of its time derivative). Actually, this relation depends on the initial values of the complex function n . 
One can suitably fix these, without loss of generality, so that & n (vo) = and 0,1(770) = — i [33|]. Then, one gets 

An = ■Z [Wnfn(Vo) + ifn(vo)}- (39) 

ZUJ n 

In the free massless case (regime which the system is expected to approach as n — > oo), one would have Qn(v) = ~ z ' 
Vti. It is convenient to reexpress 0„ so that this contribution is extracted: 

e»fa) = -*fa - vo) + f n y^-dfj. (40) 



The complex function W n that we have introduced satisfies W n (rj ) = 0. Eq. ([25]) now transforms into the Ricatti 
equation 

W n = iuj n r n -s n + (2ioj n - r n )W n - W%. (41) 

A careful but straightforward asymptotic analysis of this equation in the large-n limit, similar to that presented in 
Ref. 20], leads one to conclude that its solutions are at most of the order of n -1 . Neglecting explicitly subdoininant 
terms in the right-hand side of Eq. (|4ip , we arrive at 

W n = -s + 2iu n W n + 0{n- r ), (42) 

where s — linin-^oo s n denotes the leading contribution to the function s n , up to corrections of order 0{n~ 2 ). Note 
that s is well defined [see Eq. (|29b[) ]. With the initial condition W n (r]o) = 0, we obtain 

W n (r)) = —- s(v) - s(J7o)e 2iw " ( "- 7 ' o) -— / s(fj)e- 2iu ^ df] + 0{ n - 2 ). (43) 



2iuj r , 



2iu> 



Assuming that the function s is differentiable and its derivative s is integrable in every closed subinterval of the 
(possibly unbounded) interval of definition of the time 77, we then see that the solution is indeed of the order of n _1 . 
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On the other hand, a simple calculation shows that the coefficients a n and /?„ have the following expressions in 
terms of the functions W n and 9 n : 

a n (v,Vo)= 1+pM {[l + D n (r 1 )][l + B n ( m )] - D* n (r,)B^( m )j, (44a) 

n (V, Vo) = i+ p ^ o) { i 1 + A,fo)] 5 "fa) ~ D *n(v) [1 + B;fa)] } , (44b) 

where 

S n = J L - ^ , D n = B* n + 1(1 +p„) — , (45) 

and p n and g„ are given in Eqs. ((251) . 

Using these relations, it can be easily checked that the coefficients f3 n are of the asymptotic order of n~ 2 . As 
a result, the sequence {n(3 n } n turns out to be square summable, and the dynamics is unitarily implement able. In 
comparison with the situation discussed in Ref. [20j , the asymptotic corrections to the equations of motion provided 
by the functions p n , q n , r n , or the 0(n~~ 2 ) contributions to s n , are subdominant enough so as to preserve the unitary 
implementability of the dynamics. 

B. Equivalence of the invariant representations with unitary dynamics 

The quantum representation we have presented is not the only SO(4)-invariant one in which the dynamics is 
implemented by a unitary operator. Nonetheless, we can prove that any other representation with an SO(4)-invariant 
vacuum state and unitary dynamics must be unitarily equivalent to the one given above. 

In the following, we will call invariant the complex structures that commute with the action of the group of SO (4) 
transformations. Actually, an invariant complex structure cannot mix modes differing in any of the labels n, I, or m. 
Firstly, since the eigenspaces of the Laplace-Beltrami operator on S 3 are actually irreducible representations of the 
SO(4) group, such eigenspaces must be preserved by any invariant complex structure J. In consequence, J has the 
general form J = @ n J„, where each complex structure J„ acts exclusively on the subspace F n ®V n of the phase space. 
Furthermore, J„ consists of four components, namely, the maps J^ , J^ T ', J^^, and J^ v , each one connecting the 
spaces indicated in the superscript — for instance, J^ V connects V n with T n . As T n and V n transform with the same 
irreducible representation of SO(4), each of these four components transform with the same representation. To be 
invariant, they all must commute with all the elements of the SO (4) group. That being the case, Schur's lemma J34( 
implies that the four maps are proportional to the identity one, and hence they must respect the labels I and m as 
well [201 ■ Therefore, one can conclude that the matrix representation of J n must be block diagonal, with each 2x2 
block j„ acting just on the linear span of {f n , if/ }, with fixed n, although its explicit form depends only on n (or, 
equivalently, on the eigenvalue of the Laplace-Beltrami operator), but not on I and m. 

If we consider now a generic invariant complex structure J, we can assign annihilation and creation variables to it. 
This new basis on phase space can be related to the old one, formed by the corresponding variables of the complex 
structure J , by a symplectomorphism K, with a block diagonal matrix representation of the same type as J |16ll20j. 
We parametrize its blocks K, n in the form 

£„ = ( \* „* ) i (46) 



a: 



Z 1 > 



where we have used again that all the maps under consideration are real, and we have \n n \ 2 — |A n | 2 = 1, so that JC 
preserves the symplectic structure. Then, employing a change of basis of this kind, the relation between a generic 
invariant complex structure and the one of reference, Jo, is J — K-JqJC -1 or, block by block, j n — K. n joJC~ . 

Let us then assume that we are given an invariant complex structure J, distinct from Jo, which allows for a natural 
implementation of the dynamics. We will now show that any such complex structure gives rise to a representation 
which is unitarily equivalent to the one defined by Jq. The necessary and sufficient condition for the unitary imple- 
mentation of the symplectomorphism U with respect to this new complex structure is that the transformation K~ X WQ 
(obtained with a change of basis as discussed above) be unitarily implementable with respect to the original complex 
structure Jq |16| . Of course, the matrix representation of K~ x bllC also consists of blocks, which we can write as 



/ 



Pi 



K n U n ICn - ( /gj\* ( a J\* ) ) (47) 
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where 



J 



\ K n) Pn 



A„^„ 



2iK* n \ n ^s(a r , 



(48) 



The symbol 9 stands for the imaginary part. We will not need an explicit expression for a„, because the unitary 
implementation of U with respect to J relies only on the square summability of the sequence {nf3^} n . In the 
following, we will show that this square summability (at all times rj) implies indeed that of {n\ n } n , which is precisely 
the necessary and sufficient condition for the equivalence of the invariant complex structure J and of the reference 
one, Jo j the relation between them being given by a unitary transformation. 

If n/3^(?7,77o) forms a square-summable sequence, at least for rj in a certain time interval / around the arbitrary 
initial time i] , the same will happen with nf3n(r):Vo)/(Kn) 2 , because |k„| 2 = 1 + |A„| 2 > 1. We can subtract to each 
term of this latter sequence the contributions that we know that are square summable (i.e., any term of the order of 
n" 1 or smaller), so as to obtain another square-summable sequence. From Eq. (|44ap . it is not difficult to see that 
a n = e"" e " -f 0(n~ 2 ), and introducing the asymptotic limit of W n (|43|) in Eq. (|40|) . one concludes that the sequence 



n — sin 



U n {V - Vo) 



v s(fj) 



dq 



(49) 



'jy 2w n 

must be square summable in /. Let Z(r\) be the sum of its terms squared. If Z were integrable in /, we might write 



A 



M 

V n 2 

n— no 



M 

<E- 2 



drj sin 






HO 



2uj„ 



< j dr)Z(r)), 



(50) 



where M is any integer bigger than no, and this latter integer is chosen so that the integral over / in the second term 
of the above inequality is bounded from below by a strictly positive quantity A„ for all n > uq . This is always 
possible because, when b->oo, the last summand in the argument of the sine tends to zero and the integral can 
be seen to approach half the Lebesgue measure of the interval / [20J. Hence, the sequence {nX n / K n } n would turn 
out to be square summable, as the partial sums of its square elements form a non-decreasing sequence that would be 
bounded from above. 

Even if Z were not integrable, since it is the limit of a sequence of (partial sums of) measurable functions, and 
therefore it is measurable itself, Luzin's theorem [35j ensures that, for all 5 > 0, there exists a set 1$ C I such that 
the measure of I \ Is is less than S and the restriction Z$ of Z to Is is a continuous function. Using this fact, we 
can restrict the integrals of expression (I5TI1) to the set Is- A lower bound for the integral of the squared sine is then 
A„ — S, which is still strictly positive provided that we choose a sufficiently small 5. Therefore, the sequence formed 
by n\ n j K n turns out to be square summable in this situation as well, for similar reasons to those explained above. 

Bearing in mind that |k„| 2 — |A„| 2 = 1, it is easy to see that the square summability of {n\ n /K n } n implies that A„ 
tends to zero and n n to the unity when n approaches infinity. Consequently, the sequence formed by K n is bounded 
from above, and thus l/n n is bounded from below. It then follows straightforwardly that the sequence {n\ n } n must 
be square summable. In this way, we conclude that the symplectomorphism /C, relating two SO(4)-invariant complex 
structures, is indeed unitarily implementable. Therefore, all invariant complex structures which allow a unitary 
dynamics are unitarily equivalent. This proves that the criteria of symmetry invariance and unitary evolution that 
we have imposed select a unique equivalence class of Fock representations. Again, we see that the modifications of 
the field dynamics that appear in the equations of motion (p25")) and (|28p with respect to those analyzed in Ref. [2fJ 
do not alter the result of uniqueness of the Fock representation, essentially because the asymptotic behavior of the 
Bogoliubov coefficients is changed only at the order of -nT 2 or smaller. 



C. Uniqueness of the field description 



We have obtained a Fock quantization for the perturbation of the scalar field in terms of the modes f n and their 
conjugate momenta. This quantization is characterized by a complex structure with SO(4) symmetry and unitary 
dynamics. In addition, we have proven that any other complex structure with the same properties leads to a unitarily 
equivalent representation of the CCR's associated with our choice of a canonical pair of field variables. However, there 
is still the possibility that a different choice of field variables might result in a different quantization. Specifically, 
our field description was reached after performing a particular scaling by a time-dependent function. The question 
arises of whether, by performing a different time-dependent scaling, one might attain a distinct field description 
in which there might exist an invariant Fock representation of the CCR's leading to a unitary evolution and such 
that the corresponding quantum theory could not be related with the one that we have constructed by a unitary 
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transformation. It is worth recalling that a time-dependent scaling of the field can be completed into a canonical 
transformation by introducing the inverse scaling of the field momentum. One may also allow for a time-dependent, 
linear contribution of the field configuration to the new momentum while respecting the linearity and locality of the 
system. Since this canonical transformation depends on time, the dynamics in the new field description differs from 
the original one. Therefore, in principle, the criterion of unitary implcmcntability of the evolution imposes different 
conditions on the diverse field descriptions that can be obtained in this way. In this section, we want to prove that, 
in fact, there is no freedom to change the field description if the criteria of SO (4) invariance and unitary dynamics 
are to be fulfilled. 

Mode by mode, the class of time-dependent canonical transformation that we are considering has the form: 

f n (v) = F(v)f n (v), (51a) 

*/>) = y^fsM + G(v)Uv)- (sib) 

It is important to emphasize that the time-dependent factors in this linear transformation are the same for all the 
modes. Thus, the canonical transformation is local in terms of the configuration and momentum field variables. 
Besides, we assume that the functions F and G are real, twice differentiable, and that F does not vanish anywhere, so 
that we neither spoil the differential structure formulation of the field theory nor introduce artificial singularities. If F 
and G were constant, the quantum representation for the new field variables would be the same as for the barred ones, 
with the new canonical pair straightforwardly obtained from the original one by mere constant linear combinations. 
Incidentally, this allows us to fix the initial conditions so that the two canonical pairs coincide at the time 770, that is, 
we can set F(rj ) = 1 and G(r] ) = with no loss of generality. 

Thus, we are interested only in time-dependent transformations, which affect the dynamics. In fact, if one defines 
annihilation and creation variables combining the new variables in a way analogous to Eq. (|30p . their evolution can 
be still expressed as in Eq. ([34]) , except that now the Bogoliubov coefficients are [lj| [21j : 

a n (v,Vo) = F + (r))a n (r],r)o) + F-(r))P*(r),T) ) + - — G(T])[a n (r),Vo) + Pn(v,Vo)], (52a) 

LLO n 

$n{v,Vo) = F+(v)(3n(v,Vo) + F-(r])a* n (r],r]o) + — G(r])[p n (r],Vo) + <(l,»»)]i (52b) 

where 2F± — F ± 1/F. Remarkably, the corresponding symplectomorphism U is not unitarily implementable with 
respect to Jo in any time interval. Furthermore, it has been shown [21J that, actually, there exists no SO(4)- 
invariant complex structure in whose corresponding representation U is implemented by a unitary operator, unless 
the transformation (|5Tj) is the trivial identity. We briefly sketch the proof here. 

Let us suppose that there is such an invariant complex structure J. It will be related to Jo by a symplectomor- 
phism /C, as we have seen above. Then, if the dynamics of the new field description of the system admits a unitary 
implementation with respect to J, K,~ X UK, must be unitarily implementable with respect to Jo- Parametrizing the 
blocks K, n of the matrix representation of /C as in Eq. (|46j). direct arguments show the square summability of the 
sequence of elements 77/3^(77,770) (at all times 77 in the considered interval I), where the Bogoliubov coefficients $£ are 
given by 

& = «) 2 p n - \ 2 J* n + 2<A„3(a„). (53) 

If this is the case, the coefficients /3^/(k*) 2 must clearly tend to zero as n — > 00. Using Eqs. (1521) . together with the 
asymptotic limits of a n and /3„, and the fact that w n — n + 0(n~ l ), we arrive at the conclusion that the sequence of 
elements 

A 



e in())-))o) _ I ^li ] e ~in(v-Vo) 



F_fa)-2»-fisin[n(tj-»to)]F + fa) (54) 



K 



must tend to zero at large n. We can fix a time of the form r\ — 7/0 + 2nq/p in the interval /, with q in N and p 
in N + , so as to simplify the above expression, and consider only the subsequence formed by the terms with n = mp 
(m G N + ). With this in mind, the real and imaginary parts of the quantities ([54"]) read 



1-R 



A; 



2 

mp 



\ K mp) 
2 



r I r)o + — ) , (55a) 



^mp \ „ ( _ , 27rg\ 



\ K mp) J V P J 



F- % + — . (55b) 
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Then, both sequences of elements must approach zero as m — ► oo. And this must be so for all possible choices of the 
positive integers q and p. But this is onl y p ossible if F_(j]q + 2irq/p) — 0, since it can be proven that the other factors 
cannot tend to zero simultaneously [19L |21| . Given that the set of points r/Q + 2irq/p under consideration is dense in 
the time interval I, and recalling that the function F is continuous [with F(rjo) = 1], we are led to conclude that F(r)) 
has to be the unit function. 

Substituting F(rf) = 1 in the terms (|54l) . it becomes clear that A n /K* tends to zero as n increases (the sine function 
cannot tend to zero at all times in / [2l|). One can then deduce from the asymptotic behavior of the sequence 
{n/3^/(K*) 2 } n that the terms 

G(rj) - 4n^- sin [nfa - rjo)] e' in(v ' 7 ' o) (56) 

K n 

have also a zero limit. Considering the imaginary part of this sequence, one can show that nA„/fc* must also tend 
to zero (again because a sequence of functions of the type sin [71(77 ~~ Vo)] sm \p(jl ~ Vo) + £»] cannot tend to zero at 
all times |21j). It follows that the sequence of terms (|56|) can have a zero limit if and only if the function G vanishes 
identically. Thus, the transformation (1511) is actually the identity; only in this case can U be implemented in a unitary 
way. 

In summary, a time-dependent scaling of the field perturbation other than the Mukhanov scaling prevents the 
unitary implementation of the evolution. Note that, for instance, this is not the scaling adopted in Ref. [22]. The 
requirements of unitary quantum dynamics and SO(4) symmetry select a preferred choice of canonical pair of field 
variables and, as we have proven, a unique unitary equivalence class of Fock representations for them. 

D. An equivalent quantization 

While implementing the canonical transformation (|21|) in Sec. UH we took care that the coupling terms between 
the scalar field perturbation and its canonical momentum in the new Hamiltonian decreased sufficiently fast in the 
limit of large n. This issue was key to the unitary implementation of the dynamics in the Fock representation 
introduced in Sec. IIII Al Of course, it is possible to perform a canonical transformation that removes those coupling 
terms completely. The price to pay is that each mode of the perturbation has to be treated separately, instead of 
transforming the perturbation as a whole. In other words, one has to perform a canonical transformation that is 
defined in a different way for each mode, and therefore cannot be defined locally in terms of the original field and 
its momentum. We address now this alternate kind of transformations. It is important, however, not to spoil the 
linearity of the system, at least. For this reason, we consider linear canonical transformations of the type (|51[) . but 
with mode-dependent functions F n and G n . It is easy to see that the election 

F n = —^= (57a) 



drops the term linear in the first derivative from the equation of motion, which can be written in the form (|28|l . 
with r n — 0, whereas s n is still given by Eq. (|29b|) . up to corrections of the order of n~ 2 . In other words, the new 
configuration field variables obey equations very similar to those for the modes of a Klein-Gordon field, except that 
the mass is a time function with mode-dependent, subdominant contributions. We complete the transformation with 

1 / 1 W n \ 

G n = —== E^ - -^^- . (57b) 

Then, the new momentum variables satisfy exactly the relation -ft * — f n . In this way, we totally eliminate the 
configuration-momentum coupling terms of the field perturbations in the Hamiltonian. 

With the system in this form, the results for the case of a scalar field with time-dependent mass on S 3 |20J can be 
applied in a straightforward way. Accordingly, there exists a unique unitary equivalence class of Fock representations 
of the CCR's for the new canonical pair of field variables which possesses the desired properties of unitary dynamics 
and SO (4) symmetry. A representative of that class is the quantization that can be constructed in terms of the 
annihilation and creation variables defined by combining the new configuration and momenta variables in the same 
way as we did with the barred ones in Eq. (|30l) . A block-diagonal symplectomorphism /C relates the new annihilation 
and creation variables with the ones previously defined. Of course, its blocks JC n can again be parametrized as in 
Eq. (|46|) . and a calculation shows that 

A " = K F --£) + i Gn - (58) 
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Recalling that E™^ = 1 + 0{n~ 2 ) and that EJ n decreases asymptotically as n~ 2 , according to Eqs. (|22|) . it is 
easy to check that the sequence {n\ n } n is square summable. This is precisely the condition for K. to be unitarily 
implementable with respect to Jo- Therefore, the new quantization is unitarily equivalent to the one discussed in the 
preceding subsections. 

Thus, we obtain an equivalent quantum theory for the scalar field perturbation even though we transform each 
mode individually, a fact which prevents us from employing directly the uniqueness results explained above about 
changes of field description resulting from canonical transformations that vary on time. It is crucial to emphasize that 
the uniqueness of the quantization is granted only because we have started from the massive scalar field and we have 
considered time-dependent scalings of it, showing that there exists exclusively one which satisfies our criteria. For 
that particular field description, we have demonstrated the uniqueness of the Fock representation, and then we have 
shown that one can carry out a unitary transformation in that quantum theory so as to absorb all coupling terms 
between configuration and momentum field variables in the system. Had we started with the quantum representation 
introduced in this subsection and admitted the possibility of performing time-dependent canonical transformations 
mode by mode, in principle, we would not have been able to remove the ambiguity in the choice of a specific field 
description. 



IV. GAUGE-INVARIANT VARIABLES 

In the study of cosmological perturbations, gauge- invariant quantities are usually employed to describe the physics in 
a consistent manner, independent of the identification of the spacetime and its matter content modulo diffeomorphisms, 
and insensitive to the specification of a particular gauge. In this section, we will construct a canonical pair of such 
quantities in terms of the canonical pair of field variables used so far, and discuss their relation with the set of 
well-known gauge-invariant variables introduced by Bardeen [5| . 

Let us consider a gauge transformation of the type x^ H > x^ + e£^, where e£ M represents a small displacement that 
we will treat as a perturbation (e being the perturbative parameter introduced in Sec. UH). The covariant counterpart 
of £ M can be decomposed using hyperspherical harmonics: 

e = V2w 2 7V ^T rl Q", & = V27rtr 2 e° £ L n (PO n , ( 59 ) 

n n 

with T n and L n being functions only of time. Genuine vector harmonics do not need to be included, as they would 
only affect the genuine vector perturbations of the metric [5(, which we have ignored. 

Under the considered transformation, the background spacetime is understood to remain the same, but the per- 
turbative coefficients appearing in Eqs. (j6|) change. Including corrections up to the order of e at most, the change is: 

a n i-» a n + e~ a (aT n + ±L n ) , (60a) 



3' 



(60b) 

f n ^f n + e~ a <pT n , (60c) 

9n i-> 9n + e- a f n , (60d) 

k n ^k n - N e~ a (uj 2 n T n +L n - aL n ) . (60e) 

Recall that the dot stands for the derivative with respect to the conformal time. It is easy to combine the above 
variables to construct gauge-invariant ones at the considered perturbative order. Two independent examples of these 
are 

(61a) 

(61b) 

where we have introduced a function proportional to the energy density of the background scalar field, 

£o = ^(e- 2 V + ™V)- (62) 
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The quantities that we have defined correspond to two of the gauge-invariant scalar quantities introduced by Bardeen, 
namely, the energy density and matter velocity perturbations 5]. In the longitudinal gauge, and using the vari- 
ables (fTT)l and (|2~Tj) . their expressions reduce to 

C = e -^~ M/„ - e~ 2a (27f a7 r v - e m m 2 £) f n ] + <3(e 2 ), (63a) 



v s n = vV - 1 — h + 0(e 2 ). (63b) 

Obviously, any combination of these invariant quantities and the background variables is also gauge-invariant. Let us 
define the functions 

^ n = —^= e — E E™, (64a) 

Vn 2 - 4 7i> 



n *- = - V -5 Te~°^v< ( 27r «^ " e ba m*<p)* n . (64b) 

in -1 -k v 

These quantities satisfy the relation ty n = H^ n + 0(e), as can be seen from the equation 

fj ft 2 4 

^(e 3 "EoC) « --j==e- 3 <*ny n + 0(e), (65) 

which holds for any isotropic matter content |5j, as it happens to be our case. In terms of f n and n i , they are given 
by 

tt n = r L— (n fn + x/n) + 0(e 2 ), (66a) 

Vn — 4 

Ihfr„= JL^ (7f /ri+x / T1 )-v / » 2 -4/ Tt + 0(e 2 ), (66b) 

where 

X = - — (2^ a7V -e 6a m 2 ^). (67) 

Using these expressions, it is easy to see that ^ n and n*^ are in fact canonically conjugate on the complete phase 
space of our reduced system (including the homogeneous degrees of freedom), as the notation suggests. In addition, 
the functions $„ satisfy the same kind of dynamical equations as the modes of a scalar field with time-dependent 
mass: 

*™ + (n 2 - 4 - 3e' 4c V 2 - ± - X 2 ) *n = 0(e). (68) 

Therefore, we can construct a Fock representation with SO(4) symmetry and unitary dynamics by defining annihilation 
and creation variables in the usual way [20J: 

a * n \ _ * f^n l\ ( ^n | /gg\ 



a%„) \TS^i V^n _ V \ n * 



Any other representation with the same properties must be unitarily equivalent to this one |20| . 

The canonical transformation that leads to the variables (^> n , n* n ) is more general than the ones we have considered 
in Sec. 1111 D[ for the new coordinates have also contributions from the old momenta. Nevertheless, it does not mix 
the labels n, I, and m, so the symplectomorphism /C connecting the two bases of annihilation and creation variables is 
still block diagonal, with blocks /C„ of the form (|4"6"|) . The antilinear part of this symplectomorphism is characterized 
by the coefficients 

A, * 2 +3 f70) 

" 2v^2^7vn2— T v ' 

Therefore, the transformation can be implemented by a unitary operator in the quantum theory, as the sequence of 
elements n\ n (where we are taking into account the degeneracy of the eigenspaces of the Laplacc-Bcltrami operator) 
is certainly square summable at all times, provided tt v does not vanish and hence x is wen defined. Thus, with this 
caveat, we could alternatively quantize the perturbation in terms of the gauge- invariant variables (^„,,II^ ra ) and 
obtain a quantum theory which is unitarily equivalent to the one which is picked up by our criteria of symmetry 
invariance and unitary evolution, circumventing in this way any gauge dependence. 
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V. CONCLUSION 

With our analysis, we have shown in a specific example that the uniqueness results for the quantization of a scalar 
field with time-dependent mass in a compact spatial manifold |2(j, [2l| can be generalized to situations in which the 
dynamical equations of the system have additional corrective terms, as long as they present a suitable subdominant 
behavior in the large-n limit or, equivalently, in the asymptotic limit of large eigenvalues ui n of the Laplace-Beltrami 
operator. In particular, we have studied the quantization of the scalar perturbations of a massive field coupled to the 
gravitational field in an FLRW model. We have considered the case in which the spatial sections are homeomorphic to 
three-spheres, as a standard example of compact topology, in which infrared divergences do not appear. Nonetheless, 
the conclusions should be generalizable to other compact topologies, like e.g. the case of flat sections with the topology 
of a three-torus, in view of the results developed recently in Ref. [36{ . 

The local gauge degrees of freedom have been removed with two alternate sets of gauge-fixing conditions: in the 
main body of the article, we have adopted the longitudinal gauge, in which the shift vector vanishes and the three- 
metric of the spatial sections is conformal to that of the three-sphere, whereas in the Appendix we have imposed 
the homogeneity of the spatial sections. After these gauge fixations, the remaining local degrees of freedom can be 
identified with those of the perturbation of the matter field. The uniqueness results of Refs. [20, [21| can be extended 
to both of these gauge-fixed systems. The modifications in the second-order field equation of the Klein-Gordon type, 
in the form of a dissipative term and of a mode- dependent correction to the mass, as well as the changes in the relation 
between the field momentum and the time derivative of the field configuration, turn out to be of the order of n -2 or 
smaller. These modifications do not alter the dynamical behavior of the modes in the ultraviolet limit (and, a fortiori, 
of the Bogoliubov coefficients) up to corrections 0(n~ 2 ), which is precisely the order that can be neglected in the 
proof of the uniqueness of the Fock quantization of the system, both concerning the adoption of a field description 
and of a Fock representation for it. 

In such a quantization, the homogeneous degrees of freedom have been treated classically, thus neglecting the 
quantum effects arising from them. This is usually expected to be a good approximation away from the Planck 
regime. Since the system has been linearized, the perturbation of the scalar field admits a standard Fock quantization. 
We have restricted to complex structures that share the SO (4) symmetry of the spatial sections — hence defining an 
SO(4)-invariant vacuum state — and permit a unitary implementation of the field dynamics. These two properties 
suffice to select a unique unitary equivalence class of representations for a suitable choice of the canonical pair of field 
variables. A representative of this class is the quantization constructed with the complex structure that would be 
natural in the case of a massless field on the three-sphere. 

Furthermore, if one explores other field descriptions with distinct canonical pairs, attainable by means of a time- 
dependent scaling of the field after completing it into a linear canonical transformation, no additional Fock repre- 
sentation can be found that meets the two requirements simultaneously. The correct scaling of the field has been 
implemented in Sec. Ill Bi as a canonical transformation on the complete phase space of the system, including the ho- 
mogeneous degrees of freedom. In addition, we have determined more general, mode-by-mode time-dependent linear 
canonical transformations that allow us to reformulate the system in the genuine form of what would be a scalar field 
with time- varying mass on the three-sphere, with no subdominant corrections, and for which one can apply previous 
uniqueness results about the Fock representation of the OCR's (based also on the same criteria) [2(| . In this way, we 
have obtained two new quantum Fock descriptions of the perturbations. Actually, one of these formulations is made in 
terms of gauge-invariant quantities: the energy density and matter velocity perturbations introduced by Bardeen [5j, 
whose definition is independent of the gauge. In particular, our analysis has led us to establish a canonical structure 
for gauge-invariant variables. We have proven that the two mentioned, alternate Fock quantizations of the system are 
unitarily equivalent to the one picked up by the criteria of symmetry invariance and unitary evolution for the scalar 
perturbations of the massive field. We also notice that the canonical pair of gauge invariants that we have found 
differs from others studied in previous works [33, [38|. Besides, let us mention that it is possible to show that this 
pair can be related to the Mukhanov-Sasaki variable (for this non-flat spatial case) together with a suitable conjugate 
momentum by means of a transformation that is unitarily implement able, therefore leading to a quantization which 
is equivalent as well [39{ . 

It is worth emphasizing that it is only for the massive scalar field that we can naturally consider time-dependent 
scalings that absorb part of the time variation of the background, while respecting the locality and the linearity of the 
system. Mode-by-mode transformations are defined in general in a non-local way. Thus, starting with the massive 
field we can select a unique scaling and choose a privileged canonical pair by demanding our criteria of vacuum 
invariance and unitary dynamics. For the corresponding OCR's, the same criteria provide a unique, preferred Fock 
representation. And, from this unique Fock quantization, a unitary transformation leads to an alternate quantization, 
adapted to the mentioned canonical pair of Bardeen's gauge invariants. The consistency in the application of our 
criteria to select the quantum theory is complete, for we have further shown that the Fock quantization determined 
in this way is precisely the unique one which implements the evolution of those gauge invariants in a unitary way, 
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while preserving the S0(4) symmetry. Note also that it is the fact that we have begun with the quantization of the 
massive scalar field what allows us to choose exactly that canonical pair of Bardeen's potentials and not any other 
one. 

The robustness that these uniqueness results confer to the physical predictions of the quantization is remarkable. 
The unitary implementation of the dynamics, together with the symmetries of the spatial sections, prove to be 
powerful criteria to select not only preferred canonical field variables, but also a preferred Fock representation for 
them. If one accepts these criteria, the conclusions of the present work have implications in the study of structure 
formation and cosmological inflation, since they pick up a quantization of the (scalar) cosmological perturbations. 
One particularly interesting application is found in the so-called hybrid quantization approach to cosmology. This 
approach, which combines the Fock quantization of the inhomogeneities with the polymer quantization of the zero 
modes of the system [4(| , allows one to include local degrees of freedom in Loop Quantum Cosmology [30, 5l], |42| ■ 
The completion of this program in the case of the scalar perturbations of an FLRW spacetime will be the subject of 
future research I43I1. 
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Appendix A: A different gauge 

In this appendix, we choose gauge- fixing conditions different from those given by Eqs. (|12l) . More precisely, we 
impose 

a n = = b n , (Al) 

so that the perturbation of the three-metric vanishes. In order to check the validity of these conditions, we calculate 
their Poisson brackets with the constraints that are going to be fixed, 

A Apt ({ a n,H n i} {b n ,H2}\_ 1 4 



{a ni H™} {b n ,H£}J -3 

which becomes zero exclusively when ir a = 0. As far as that point is eluded (or, rather, that section on phase space), 
the fixation of the gauge is admissible. We still have to require its stability under the dynamics, which implies 

= {a ni H} « -N e~ 3a (n an +g n Tr a ) - ^e~ a k n + 0(e), (A3a) 

n 2 — 1 1 

= {b n ,H} » Aoe" 3 "-^— -n bn + -e~ a k n + 0(e). (A3b) 

The latter equations can be solved for the Lagrange multipliers g n and k n . In turn, the values of 7r ari and 7Tf, w can be 
obtained from the constraints H™ and 7?,", together with the conditions (|A1|) . One gets 

n an = — {~K V -Kf n + e 6a m 2 ipf n ) , TT bn = — [7i>7r/ n - (371-^ - e 6a m 2 ip) f n ] . (A4) 

7T a 7T a 

After reduction of the system, the Hamiltonian has the same structure as in Eq. (TT8"|) . The zeroth-order Hamiltonian 
is still given by Eq. ([8]), while the second-order Hamiltonian can be written in the form (fl~9f with the coefficients 






2 



En - = e ~ Ja V + —±i V (A5a) 



E}\ = -3e- 2a 



7r„ n z - 4 -ni v ' 



(A5b) 



E ff = e (n -1 



-2a 



6a 2, 



9< - 6e Da mV 



+ e 4Q m 2 + 



3e- 2a 1 /„ 



6a 2 V 

em y>J 
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(A5c) 



Let us now perform the Mukhanov scaling of the field, completing it into the following canonical transformation 
(up to the considered perturbative order): 



q = q — 



e- 1 £& 






./ n e In 



^ 2 E 



- /n7T/„ 



T/n 



■ 7I"a /, 



3-^+TTa /, 



(A6a) 
(A6b) 



whereas tt^, is left unchanged. This transformation does not alter the structure of the Hamiltonian, but absorbs the 
dominant contribution in the large-n limit from the terms that couple the matter-field configuration modes with their 
momenta. The new coefficients of the second-order Hamiltonian are 



El 



4*?' 



m, = ^^h< 



j f* 



Ani 



2-KaTTt, 



e m m 2 Cp 



(A7a) 
(A7b) 



E 



if 



- + e 2a m 2 



37T 

2tt 



— e 

2 



-4a 



(*% - 30^ 2 ) 



(n*-4) ffl 



27r a 7r v + e oa m^ip 



— e 
2 



-4a 



27^f + 3e 6a m 2 £ 

TTa 



^-£3^-1 



(A7c) 



Applying Hamilton's equations in conformal time, one finds that the momentum ttf is still given by Eq. (|25|) . 
except that now 



37T 2 



Pi, 



-3e 



-2a„ 



3(tt3 /7f g ) - 27r 57 r v + e 6a m 2 £ 



(n 2 - 4)^ + 3^' - - (n 2_ 4) ^2 +37r 2 

As for the equation of motion of f n , Eq. (|28|) is still valid as long as we redefine the coefficients in it to be 



(A8) 



_ o -2a *V 37F " ~ 3# g^ + e 6a ^V(27Ta ~ 3^) + e 4g 7T ¥ 



(n2-4)7f 2 +37r2 



s„ = - + e 2a m 2 



2tt1 2 



30tt£ + 27^ 



2 C 



-4a 



24e w Wu-£ - 3e oa m a w a 3^-1 



(A9a) 

+ 0(n- 2 ). 
(A9b) 



Note that the expressions for p„, g„, r ni and s n in the new gauge have the same asymptotic behavior as in the gauge 
adopted in the main body of the article. Owing to this fact, it is not difficult to go through the same arguments and 
check that the results of Sec. IIIII apply as well in this case. 

The quantities defined via Eqs. ([M]) continue to satisfy a relation of the form ^ n = H^ n at leading perturbative 
order, because Eq. (|65[) is expressed in terms of gauge-invariant quantities. Again, they are canonically conjugate, 
and the equations of motion for ty n are analogous to those for the modes of a Klein-Gordon field with time-dependent 
mass. Hence, the corresponding OCR's admit a Fock representation which implements the dynamics in a unitary 
manner. For this, we can introduce two sets of annihilation and creation variables, as in Eqs. (|30|) and (|69|) . The 
symplectomorphism that relates these variables with the ones adopted in a quantization of the kind discussed in the 
main text is block diagonal and can be checked to be unitarily implementable, since the coefficients A„ of its antilincar 
part still have the form (|70p . but with the function \ replaced with 



X 



6a™2 ~ 



3^ - 27f a 7iv + e^m'ip 



(A10) 



This new function is well defined if so is the selected gauge (so that, at leading order in the perturbative expansion, 
the vanishing of 7f g is obviated) and provided, once more, that ir v differs from zero. 
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In conclusion, we see that the results of Sec. IIVI are reproduced in this gauge, a fact that proves their robustness. 
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